No evidence for spontaneous orbital currents in finite size studies 
of three-band models for CuO planes 
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We have numerically evaluated the current-current correlations for three-band models of the CuO 
planes in high-T c superconductors at hole doping x = 1/8. The results show no evidence for 
the orbital current patterns proposed by Varma. If such patterns exist, the associated energy is 
estimated to be smaller than 5 meV per link even if e p — £d = 0. Assuming that the three-band 
models are adequate, quantum critical fluctuations of such patterns hence cannot be responsible for 
phenomena occurring at significantly higher energies, such as superconductivity or the anomalous 
properties of the pseudogap phase. 

PACS numbers: 74.20.Mn, 74.72.-h, 74.20.-z 



The problem of high-T c superconductivity has been 
something like a holy grail to the field of condensed mat- 
ter physics for the past two decades [l[. In the words 
of R.B. Laughlin, it has not been a fight, but a war. 
It has been a traumatic experience for some of those 
involved, but has also led to a plethora of new devel- 
opments extending far beyond the field. Many of the 
experimental techniques used to study the systems, like 
angle resolved photo emission spectroscopy (ARPES) or 
scanning tunneling microscopy (STM), have undergone 
revolutions with regard to resolution and data process- 
ing. The theory of superconductivity in high-T c cuprates 
has been found as many times as victories must have 
been proclaimed in civil wars, but while individuals be- 
lieve to have the theory, there is no consensus what the 
theory should be. Many ideas, even though too general 
to qualify as complete theories of the cuprates, have in- 
spired a vast amount of research in both high-T c and 
other areas. Most prominently among them are the no- 
tions of a resonating valence bond (RVB) state the 
gauge theories of antiferromagnetism 0], and the notion 
of quantum criticality There have been, however, a 
few concise proposals which make falsifiable predictions. 
Intellectual masterpieces among them have been the the- 
ory of anyon superconductivity [H, the proposal of ki- 
netic energy savings through interlayer tunneling Q , the 
SO (5) theory of a common order parameter for supercon- 
ductivity and magnetism Q , and a more recent proposal 
that the anomalous properties of the cuprates may be 
due to quantum critical fluctuations of current patterns 
formed spontaneously in the CuO planes [§, 9]. This last 
proposal is further investigated in this Letter. 

The idea of a spontaneous symmetry breaking through 
orbital currents was, as with so many major advances 
in physics, motivated by experiment. The normal state 
of the cuprates at optimal doping shows a behavior 
which can be classified as quantum critical, and has been 
rather adequately described by a phenomenological the- 
ory called marginal Fermi liquid [id ]. This phenomenol- 
ogy suggests a quantum critical point (QCP) at a hole 



doping level of x c « 0.19, an assumption consistent with 
a significant body of experimental data [ll|, 12, 13|, 14, 
Hit ]. Critical fluctuations around this point would then 
be responsible for the anomalous properties of the pseu- 
dogap phase, and provide the pairing force responsible 
for the superconducting phase which hides the QCP. 

Interpreting the phase diagram in these terms, one is 
immediately led to ask what the phase to the left of the 
QCP, i.e., for x < x c , might be. The theory would 
require a spontaneously broken symmetry beyond the 
global U(l) symmetry broken through superconductiv- 
ity (which is often erroneously refered to as a broken 



gauge symmetry [ly]). In addition, as the fluctuations 



are assumed to determine the phase diagram up to tem- 
peratures of several hundred Kelvin, the characteristic 
energy scale of the correlations inducing this symmetry 
violation must be at least of the same order of magni- 
tude. No definitive evidence of such a broken symmetry 
has been found up to now, even though several possi- 
bilities have been suggested. These include stripes fl7| . 
a ci-density wave [lijh an d most recently a checkerboard 
charge density wave [191 ] . 

The general consensus is that the low energy sec- 
tor of the three-band Hubbard model proposed for the 
CuO planes (see (JJ) below) [20[ reduces to a one-band 
t-t'-J model, with parameters t « 0.44, t' « —0.06, 
and J « 0.128 (energies throughout this article are in 
eV) [El, 0, M, M, UM- For the undoped CuO planes, 
the formal valances are Cu 2+ and 2 ~. As the electron 
configuration of Cu atoms is [Ar] 3d 10 As 1 , this implies 
one hole per unit cell, which will predominantly occupy 
the 3d x 2_ y 2 orbital. As the onsite potential e p in the O 
2p x and 2p y orbitals relative to the Cu 3d x 2_ y 2 orbital 
is generally assumed to be of the order of e p = 3.6 (with 
£d = 0), and hence smaller than the onsite Coulomb re- 
pulsion Ud ~ 10.5 for a second hole in the 3d x 2_ y 2 orbital, 
it is clear that additional holes doped into the planes will 
primarily reside on the Oxygens. The maximal gain in 
hybridization energy is achieved by placing the additional 
hole in a combination of the surrounding O 2p x and 2p y 
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We hence undertook to estimate these couplings 
through numerical studies of finite clusters containing 8 
unit cells, i.e., 8 Cu and 16 O sites, and periodic bound- 
ary conditions (which do not frustrate but should en- 
hance the correlations). The total number of holes on 
our cluster was taken N = 9 (5 up-spins and 4 down- 
spins), corresponding to a hole doping of x = 1/8. We 
had hoped that the energy associated with a domain wall, 
which may be implemented through a twist in the bound- 
ary conditions, would provide information regarding the 
coupling aligning the orbital currents in neighboring pla- 
quets, while the splitting between the lowest energies for 
a finite system would provide an estimate for the trans- 
verse field. The result of our endeavors, however, is a 
daunting disappointment: The coupling is zero within 
the error bars of our numerical experiments. 

Let us now report our numerical studies in detail. To 
begin with, we wish to study the three-band Hubbard 
Hamiltonian H = H t + H\j with 

H t =^e p n? .-i pd ^ [4,aPj,a +p},o- d i 

i,a (i,j),<? 
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FIG. 1: Orbital current pattern proposed by Varma. 

orbitals with the same symmetry as the original hole in 
the Cu 3d x 2_ y 2 orbital, which requires antisymmetry of 
the wave function in spin space, i.e., the two holes must 
form a singlet. This picture is strongly supported by data 
from NMR [26| and even more directly from spin-resolved 
photoemission [13]. In the effective one-band t-J model 
description of the CuO planes, these singlets constitute 
the "holes" moving in a background of spin 1/2 particles 
localized at the Cu sites. 

In contrast to this picture, Varma [H, Q has proposed 
that the additional holes doped in the CuO planes do 
not hybridize into Zhang-Rice singlets, but give rise to 
circular currents on O-Cu-0 triangles, which align into 
a planar pattern as shown in Fig. [TJ He assumes that 
the inter-atomic Coulomb potential V p d is larger than 
both the hopping t p d and the onsite potential e p of the 
O 2p orbitals relative to the Cu 3d x 2_ y i orbitals, an as- jj _ jj y^ n P n P 4- jj a 
sumption which is not consistent with the values gen- P : : 

erally agreed on (see the list below fl}). Making addi- 
tional assumptions, Varma has shown that the circular 
current patterns are stabilized in a mean field solution of 
the three-band Hubbard model. The orbital current pat- 
terns break time-reversal symmetry breaking (T) and the 
discrete four-fold rotation symmetry on the lattice, but 
leave translational symmetry intact. The current pat- 
tern is assumed to disappear at a doping level of about 
x c f=a 0.19. The phenomenology of CuO superconductors, 
including the pseudogap and the marginal Fermi liquid 
phase, are assumed to result from critical fluctuations 
around this QCP, as outlined above. 

Motivated by this proposal, several experimental 
groups have looked for signatures of orbital currents or 
T violation in CuO superconductors. While there is no 
agreement between different groups regarding the mani- 
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where ( , } indicates that the sums extend over pairs of 
nearest neighbors, while d; :(T and p JiCT annihilate holes in 
Cu 3d x 2_ y 2 or O 2p orbitals, respectively Hybertsen et 
al. [23| assumed t p d = 1.5, t pp = 0.65, Ud = 10.5, U p = 4, 
V pc \ = 1.2, and e p = 3.6. 

In order to be able to diagonalize (fTJ) for 24 sites, we 
need to truncate the Hilbcrt space. A first step is to 
eliminate doubly occupied sites, which yields the effective 
three-band t-J Hamiltonian 



H eS = P G H t P G + Hj with 



(2) 



b 3 = J Pd £ (s?.s« -\)+ j pp J2 (sis* - \ 
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festation of T violation in ARPES studies 28|, 2j| , a re- 
cent neutron scattering experiment by Fauque et al. (3pl | 
indicates magnetic order within the unit cells of the CuO 
planes. Their results are consistent with Varma's pro- 
posal, and call the validity of the one-band models into 
question. 

In a recent article, Aji and Varma [3~0 | have mapped 
the four possible directions of the current patterns in each 
unit cell onto two Ising spins, and investigated the crit- 
ical fluctuations. Within this framework, the coupling 
between and the transverse fields for these Ising spins 
decide whether or under which circumstances the model 
displays long-range order in the orbital currents. 
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and the sums in H j are limited to pairs where both neigh- 
bors are occupied by holes. If Pg only eliminates config- 
urations with more than one hole on a site, the dimension 
of the St ot = \ sub-sector is with 164,745,504 still beyond 
our capabilities. We have hence implemented two further 
ways of truncating the Hubert space: (a) We limit the 
maximal number of holes allowed in the O orbitals to 
iV™ ax . (b) We limit the maximal number of CuO links 
occupied with 2 holes to N™^. For the values proposed 
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FIG. 2: Current-current correlations (jk,k+xjl,m) multiplied 
by 10 2 for the ground state of © with e p = 3.6 on a 24 
site cluster (8 Cu = open circles, 16 O = filled circles) with 
PBCs. The reference link is indicated in the top and (due to 
the PBCs) bottom left corner. 

by Hybertsen et al, either truncation should hardly af- 
fect the low energy sector. In Table U we compare the 
ground state energies and the current-current correlation 
on diagonally and maximally separated links on the torus 
for different values of N™ 3 * and iVjPjf. This compari- 
son gives us confidence that it is not unreasonable to set 
iV[™| x = 4, which reduces the Hilbcrt space dimension 
to 93,595,824. Since the truncation (b) is predominantly 
projecting out states with high kinetic energy, we expect 
that the value of e p will not affect the validity of the 
approximation. Excluding configurations with too many 
CuO links occupied by two holes should, in any occasion, 
not weaken the tendency to form orbital current patterns. 
Truncation (a), by contrast, is no longer reasonable if e p 
is small. 

Let us now turn to our results for the current-current 
correlations in the ground state (situated at the M point 
in the Brillouin zone) for our 24 site cluster with 9 holes 
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TABLE I: Ground state energies per unit cell and a current- 
current correlation for various truncations of the Hilbert 
space. 




FIG. 3: as in Fig. [2j but with e p = 0. Except for the vertical 
links, positive correlations indicate alignment with the current 
pattern shown in Fig. [TJ 

and parameters as assumed by Hybertsen, where no or- 
bital current pattern are expected. With the current op- 
erator for, e.g. an 0-0 link given by 

jk,l = H PP [plaPk,* ~ Pk,aPl,a) > ( 3 ) 

the correlations function (jk,k+xjl,m) with an 0-0 link 
as reference link is depicted in Fig. [2] As expected, the 
correlations fall off rapidly, and there is no indication of 
order. 

The crucial result is that the correlations change con- 
tinuously and only quantitatively, but not really quali- 
tatively, as e p is lowered from 3.6 to 1.8, 0.9, 0.4, and 
finally to e p = 0, where the result is shown in Fig. [3] 
(the ground state is now situated at the L point). This 
is the situation for which Varma has proposed that the 
current pattern sketched in Fig. [1] would occur. Fig. [3l 
by contrast, shows no alignment of the currents. 

The numerical experiments for the finite cluster can, 
as a matter of principle, never rule out that a symmetry 
is violated. We can use them, however, to put an up- 
per bound on the size of the spontaneous currents, and 
hence the energy associated with these currents. If a 
current pattern as sketched in Fig. [1] were to exist, the 
current correlation (jk,k+xjl,l+x) f°r links far away from 
each other in a rotationally invariant ground state should 
approach | (x\jk,k+x\x) 2 , where \x) denotes a state with 
a spontaneous current pointing in x direction (the fac- 
tor | arises because by choosing our reference link in re- 
direction, we effectively project onto two of the four pos- 
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sible directions for the current pattern) . From the values 
of 10 2 ■ (jk.k+xji,i+x) for the four horizontally connected 
links in the center of Fig. El -0.1832, +0.0158, -0.0676, 
and +0.0158, which should all be positive if a current 
pattern were present, we estimate 10 2 • (jkM+x.) 2 < 0.20 
and hence (x \jk,k+x \%) < 4 • 10 -3 as an upper bound for 
a current pattern we are unable to detect through the 
error bars of our numerical experiment. We now denote 
(x\jk,k+x\x) by j pp . 

We roughly estimate the kinetic energy e pp per link 
associated with a spontaneous current j pp of this mag- 
nitude using jpp = n p v and e pp = ^n p mv 2 with m = 
l/2f pp , where n p is the hole density on the Oxygen sites 
(?i p = 0.30 for the state of Fig. [3]), and obtain 



We further acknowledge the support of the computing 
facilities of the INT at the Forschungszcntrum Karlsruhe. 



Jpp 



4f n 



< 5 ■ 10" 



A similar analysis for the CuO links, using data not 
shown here, yields with 10 2 - (jk,k+x+y) 2 < 1.0 and hence 
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< 10 2 (there is no factor | in this case) an estimate 
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We conclude that while we cannot rule out that orbital 
current patterns exist, we can rule out that they are re- 
sponsible for the superconductivity, the properties of the 
pseudogap phase, or the anomalous normal state prop- 
erties extending up to temperatures of several hundred 
Kelvin, as the energy associated with the spontaneous 
loop currents is less than 5 meV per link if such currents 
exist. We have assumed that the CuO planes are ade- 
quately described by the three-band Hubbard model |T|) , 
but we have allowed e p to be much smaller than gen- 
erally agreed upon, and based our estimates on the ex- 
treme and to our purposes most unfavorable value e p = 0. 
(Note that the ordered antiferromagnetic phase in un- 
doped cuprates requires a finite e p .) Numerical data not 
presented here show that our conclusions remain intact 
if we set J pc j = J pp = or/and double the value of 
the repulsion Vpd, which generates the orbital currents in 
Varma's mean-field calculation. They also hold for other 
low energy states for the finite system (e.g. as situated 
at the M point in the Brillouin zone for e p = 0). 

Nonetheless, we should keep in mind that any analy- 
sis of a model can only reach a conclusion valid for this 
model. The question of whether current patterns exist 
in CuO superconductors can ultimately only be settled 
by experiment. We consider it likely, however, that an 
eventual consensus among experiments will confirm our 
conclusion. 

We wish like to thank CM. Varma, V. Aji, F. Evers 
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arship from the Studienstiftung des deutschen Volkcs. 
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